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; Abstract 

o\ ; 

In this work we present cosmological quantum solutions for all Bianchi Class 
^ . A cosmological models obtained by means of supersymmetric quantum me- 

Q | chanics . We are able to write one general expression for all bosonic com- 

ponents occuring in the Grassmann expansion of the wave function of the 
Universe for this class of models. These solutions are obtained by means of a 
more general ansatz for the so-called master equations. 
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00 1 

-st ! I. INTRODUCTION 

o ■ 

7-h ' In recent years, many particular solutions of the Wheeler-DeWitt (WDW) equation in 

standard quantum cosmology have been found for Bianchi Class A models [[j]. In general, the 
wave functions of these models have the form \l/(q M ) = W(q M ) e ±4> ^ qM - ) , where $ is a solution 
to the corresponding classical Hamilton- Jacobi equation. One can find similar solutions for 
5^ \ a WDW equation that has been derived in the bosonic sector of the heterotic string M. 

Some time ago, Graham |||4[] considered the WDW equation in the ADM formulation of the 
Bianchi type IX model and found the appropriate factorization in terms of the operators 
^ which may be considered the square roots of this equation. Graham was the first to provide 

a master equation for the Bianchi type IX models in Witten's supersymmetric quantum 
mechanics for the so-called auxiliary function f±(q M ). However, he studied only the trivial 
case f = const. Next, Obregon et al, obtained ans solved the master equation for the 
Bianchi type II models. For an application of the technique to the 2D Fokker-Planck equation 
see 0. 

Our main goal in this work, is to find the general form for all bosonic functions that 
appear in the Grassmann expansion of the wave function of the Universe, in the approach of 
Witten's supersymmetric quantum mechanics for all Bianchi Class A cosmological models. 
In the Grassmannian calculus, one introduces fermionic variable without a direct physical 
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meaning. Physical quantities of interest, such as require integration over the Grassman- 
nian variables, obtaining in this way information about their influence on the unnormalized 
probability distributions. 

We begin with metrics of the form 

ds 2 = (N 2 - N j Nj)dt 2 - 2Ni dt &j - e 2a(t) e 2/3(t) yu/ u\ (1) 

where a(t) is a scalar and /3y is a 3 x 3 matrix, and N(t) and Ni(t) are the lapse and shift 
functions, respectively. The one-forms uo l have specific expression for any particular Bianchi 
Universe ||. 

The Arnowitt-Deser-Misner (ADM) action has the form 

1 = y"(P x dx + P y dy + P Q da-Nft)dt, (2) 

where 

H = e- 3Q (-P 2 + P 2 + Pj + e 4Q V(x, y)) . (3) 

The potential term e 4ct V(x, y) = U(a,x, y) is specific of the cosmological model under con- 
sideration, being related to the function as follows |] (g M = (a,x, y)). 

U(q") = rT^r- IT-- ( 4 ) 

This turns into the Einstein-Hamilton- Jacobi equation when we replace the momentum 
P M -»■ d<$>/d<f in (D, rT = diag(-l, 1, 1). 

In this work we drop the factor e~ 3a . This corresponds to fixing the factor ordering 
before factorizing the WDW equation. The latter is got by replacing P M — > — i<9 qM in (Q). 

The work is organized as follows. In Sec. II we give the factorization of the WDW 
equation in terms of the supercharges operators and obtain the so-called master equations 
for an auxiliar function i± in general way for all Bianchi Class A cosmological models. In Sec. 
Ill we solve this master equation and present all bosonic components of the wave function 
for the same models. In Sec. IV, we write a general form for an unnormalized probability 
density, by imposing some soft of boundary conditions on the bosonic components. Sec. V 
is devoted to conclusions 



II. FACTORIZATION OF THE WDW EQUATION 

The WDW equation can be obtained as the bosonic sector of a super-Hamiltonian in 
superspace H ss by employing supercharge operators. In the 3D case, the supercharges read 



Q = V 



(5) 



and 



.9$ 
l dq u 



(6) 
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with the following algebra for the variables ip^ and ip u > 

{r, r} = v" v , {v, v] = o, {r, p} = o. (?) 

Using the representations p = 7] fiu d/d6 l/ and ip u — 0", one finds the superspace Hamiltonian 
in the form 

h. = 5 {Q,Q} = «. + §j^ [#*,*•], (8) 

where 7^ = □ + ^ (<? M ) is the standard WDW equation, □ is the 3D d'Alambertian in 
the q M coordinates with signature (- + +),{,} means the anti-commutator, and [, ], the 
commutator. Notice that the quantized super-Hamiltonian differs from the classical one H.q 
by a spin term, which vanishes in the classical limit. 

The supercharges Q, Q and the super-Hmiltonian satisfy the following algebra 

^{Q,Q} = H SS , {H SSl Q} = } {H SSl Q}=0. (9) 

Physical states are selected by the constraints 

q^ = o, g* = o, (io) 

where the following 3D Grassmann representation of the \I/ is used 

* = a + + B v e v + ^ A c x e» q v + A- e° e 1 e\ (n) 

(jl, u, A running over 0, 1, 2. 
The following ansatz 

= (12) 

in ([10|) and (0), leads to the master equation for the auxiliary function f + 

□f + -2^^^i = 0. (13) 

In addition, it is possible to show that \e^ v \ C x 9 a 9^ 9 U = C a 9°9 1 9 2 and employing the ansatz 

<9f 

C a = v a ^e* (14) 
we obtain the second master equation of the form 

□F_ + 2»T-?-S JL = 0. (15) 
Thus, (|T3| ) and ( ^5|) can be rewritten as 

The solutions for the other functions are 

A ± = q ± e^ (17) 

where q± are integration constants. 

In the following, we give the solution of fll6|) , using as toy models the Bianchi Class A 
cosmological models (l|] 
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III. QUANTUM SOLUTIONS FOR BIANCHI CLASS A MODELS 



To solve ( fPo| ) it is necessary to know the potentials U(q M ) (or $(g M )) for the Bianchi 
Class A cosmological models. These are given in table I [1 



Bianchi type 


Potential U 


$ 


I 








II 


i e m 


6 


VI h= _x 


3 C 


±1 


2(ft-ft)]e*+* 


VII h=0 


1 
3 


e 4fh _|_ e 4/3 2 _ 2e 2(/3i+/3 2 ) 


±1 


e 2 ^ 1 + e 2 ^ 2 


VIII 


1 
3 


e 4/3i + e 4/3 2 + e 4/3 3 _ 2 e 2(/3i+/3 2 ) + 2 e 2 Wi+M + 2e 2(/32+/33) 


±1 


e 2/3i _|_ e 2/3 2 _ e 2/3 3 


IX 


1 
3 


e 4/3i _|_ e 4/3 2 + e 4/3 3 + 2e 2(/3i+/3 2 ) + 2e 2(/?1+ftl) + 2e 2(/32+/33) 


±1 


e 2/?l + e 2/3 2 + e 2/? 3 



Table I 

Potential U and superpotential $ for all Bianchi Class A models, where 
Pi = a + x + y/3 y, /3 2 = a + x — V3y and /3 3 = a — 2x. 

Once f± are obtained, all the bosonic components that appear in the Grassmann expan- 
sion of the wave function ([n]) can be determined as follows. The ansatz 



f±(q M ) =W±(q' 



turns ( |I~6|) into 



whose solutions are 



(V$) 2 T n$ 



W± 



w±(q* 



# t e m " t "e =F *, 



(18) 



(19) 



(20) 



where f3± are integration constants, m M = (m 1 ,m 2 ,m 3 ) is a vector of null measure (i.e. 
— m 2 + m 2 + m 2 = 0) and t M = (—a, x, y). 
Thus, the solutions for f ± are given by 



f± = /3±e m " . (21) 
In the Table II, we write the corresponding m M t M terms for all Bianchi Class A models. 



Bianchi type 


m M P" = -m 1± a + m 2± x + m 3± y 




I 


-m 1± a + m 2± x + m 3± y 


(m 1± ,m 2± ,m 3± ) 


II 


b±(2a - x+ v^y) 


b±(-2,-l, v / 3) 


VI h= _i 


b± (a + x) 


b±(-l,l,0) 


VII h=0 


b± (a + x) 


b±(-l,l,0) 


VIII 





(0,0,0) 


IX 





(0,0,0) 



Table II 

m^t^ and terms for all Bianchi Class A models. 
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Moreover, the solutions for all bosonic components of the wave function are 



A± 


= q±e T 






Bo 


= -mi + f+e , 


C°~- 


= mi_ t_ e 


Bi 


= m 2+ f+e - *, 


C l = 


m 2 _ f_ e*, 


B 2 


= m 3 f+e"*, 


C 2 = 


m 3 _ f_e*, 



(22) 



where q±andb± are arbitrary constants and m i± are constants given in table II. 

IV. THE UNNORMALIZED PROBABILITY DENSITY |^| 2 

In this section we follow ref. to get |\l/| 2 for the wave function given by (|TT|). Integrating 
over the Grassmann variables 6^ ||, we are led to the following result 

|^| 2 = A\ A+ + B* Bo + B\ B l + B* 2 B 2 + C°* C° + C 1 * C 1 + C 2 * C 2 + A*_ A- . (23) 

Assuming that the constants appearing in (E3) are real, we obtain 



I* 



2 



q + + (mL + mi + mi) f, + q 2 _ + (mf_ + m 2 _ + m 2 J f! ^. (24) 



For example, for the Bianchi type II models, we obtain from (p4|) 

|* n | 2 = U + 8b + i%) e- 2 * + (q 2 + 8b 2 f 2 ) e 2 *, (25) 



which is Eq. (4.2) in [§. 

Now, by imposing the some sort of boundary conditions, only one kind of terms will 
remain in the wave function || . There are many possibilities, one is to demand that ^ does 
not diverge for |x|, |y| — > oo at fixed a. In such a case only the terms e - * will remain in 



Figures 1 and 2, display the behavior of the probabilitity density \^\ for the Bianchi 
type VI/^-i and VIIh=o when the parameter b+ = —1, respectively. The |\1/| 2 for the same 
models when b+ = +1 are plotted in Figures 3 and 4. The 2 behavior for the Bianchi 
type II and IX, has been given in and 0, respectively. 

V. CONCLUSION 

Our main goal in this work, was to obtain a general expression for all bosonic functions 
that appear in the expansion of the wave function of the Universe using the approach of 
Witten's supersymmetric quantum mechanics. This has been done for all Bianchi Class A 
cosmological models. In addition, we found the general form for the unnormalized probabil- 
itity density, (|24]), for any Bianchi Class A cosmological model. When we apply boundary 
conditions on the wave function, the unnormalized probability density is reduced to 



q^ + tm 2 +m 2 + m 2 ) f 2 e" 2 *. (26) 



At this step, the main contribution to l^l 2 will depend on the values of the parameter 
b+ that occurs in the function f+ and on the range used for the coordinates a and (x, y). 



5 



VI. ACKNOWLEDGMENTS 

We would like to thank H. Rosu for valuable hints. This research was supported by 
CONACyT. ERM was supported by a CONACYT student fellowship. 



6 



REFERENCES 



[1] O. Obregon and J. Socorro, Int. J. Theor. Phys. 35, 1381 (1996). For review, see A. 

Macfas, Gen. Rel. and Grav. 31, 653 (1999). 
[2] J.E. Lidsey, Phys. Rev. D 49 R599 (1994). 
[3] R. Graham, Phys. Rev. Lett. 67, 1381 (1991). 
[4] R. Graham, Phys. Rev. D 48, 1602 (1993). 

[5] O. Obregon, J. Socorro and J. Bemtez, Phys. Rev. D 47, 4471 (1993). 

[6] M.P. Ryan and L.C. Shepley, Homogeneous Relativistic Cosmologies, (Princeton Univer- 
sity Press, Princeton, NJ 1975). 

[7] H.C. Rosu, J. Socorro and O. Obregon, J. Phys. A: Math. Gen. 29, 1825 (1996). 

[8] L.D. Faddeev and A. A. Slavnov, Gauge Fields: An Introduction to Quantum Theory 
(Addison- Wesley, Reading, MA. 1991), sec. 2.5. 

[9] S.W. Hawking and D.N. Page, Phys. Rev. D 42, 2655 (1990). 



7 



Fig. 1 

At the left we show |\I/vi| 2 = q+ e~ 2 *, for the first pure bosonic term, when q + = 1, a = —6. At 
the right, the main contribution to the unnormalized probability density due to the B„ terms, 

when (3 + = 1, b + = —1, a = —6. 




Fig. 2 

At the left we show l^vnl 2 = q+ e -2 *, for the first pure bosonic term when q + = 1, a = —6. At 
the right, the main contribution to the unnormalized probability density due to the terms, 

when P + = 1, b + = —1, a = —6. 
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Fig. 3 

At the left we show l^vil 2 = e -2 *, the main contribution to the unnormalized probability 
density due to the first pure bosonic term, when q + = f , a = — 6. At the right, the corresponding 
contribution of the £>„ terms, for f3 + = 1, b + = 1, a = —6. 




Fig. 4 

At the left we show l^vnl 2 = q+e -2 *, the main contribution to the unnormalized probability 
density due to the first pure bosonic term, when q+ = l,a = —6. At the right, the corresponding 
contribution of the terms, for /?+ = 1, b + = 1, a = —6. 
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